Nonlinear optics of surface plasmon toy black holes. 
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Nonlinear optical behavior of recently introduced surface plasmon toy black hole models has been 
considered. Physical realization of these models involves droplets of liquid on a metal surface, which 
supports propagation of surface plasmons. Such droplets are shown to exhibit strongly enhanced 
nonlinear optical behavior in the frequency range near the surface plasmon resonance of a metal- 
liquid interface. This enhancement may be responsible for the missing orders of magnitude of field 
enhancement in the surface enhanced Raman scattering effect. In addition, experimental observation 
of the recently predicted optical second harmonic generation near the toy event horizon is reported. 
Finally, the possibility of toy Hawking radiation observation is discussed, and the expression for the 
effective Hawking temperature of a toy surface plasmon black hole is derived. The effective Hawking 
temperature appears to be in the range of a few hundreds Kelvin for the realistic droplet shapes, 
which may lead to a surprising conclusion that the shapes of the liquid droplets on metal substrates 
are to some degree determined by the Hawking radiation. 



PACS no.: 78.67.-n, 04.70. Bw 

The realization that solid-state toy models may help 
in an understanding of electromagnetic phenomena in 
curved space-time has led to considerable recent effort 
in developing toy models of electromagnetic [1] and sonic 
[2] black holes. In the case of media electrodynamics 
this is possible because of an analogy between the prop- 
agation of light in matter and in curved space-times: it 
is well known that the Maxwell equations in a general 
curved space-time background gik{x,t) are equivalent to 
the phenomenological Maxwell equations in the presence 
of a matter background with nontrivial electric and mag- 
netic permeability tensors €ij{x,t) and iiij{x,t) [3]. In 
this analogy, the event horizon corresponds to a surface of 
singular electric and magnetic permeabilities, so that the 
speed of light goes to zero, and light is "frozen" near such 
a surface. In the absence of established quantum gravita- 
tion theory the toy models are helpful in understanding 
electromagnetic phenomena in curved space-times, such 
as Hawking radiation [4] and the Unruh effect [5] . Unfor- 
tunately, until recently all the suggested electromagnetic 
black hole toy-models were very difficult to realize and 
study experimentally, so virtually no experimental work 
was done in this field. 

Very recently a simple surface plasmon toy black hole 
model has been suggested and observed in the experi- 
ment [6]. In addition, a surface plasmon toy model of 
a rotating black hole (Kerr metric [7]) has been intro- 
duced [8]. These models are based on the surface plas- 
mon resonance at the edge of a liquid droplet deposited 
on the metal surface, which supports propagation of sur- 
face plasmons. Surface plasmons are collective excita- 
tions of the conductive electrons and the electromag- 
netic field [9]. They exist in "curved three-dimensional 
space-times" defined by the shape of the metal-dielectric 
interface. Since in many experimental geometries sur- 
face plasmons are weakly coupled to the outside world 



(to free-space photons, phonons, single-electron excita- 
tions, etc.) it is reasonable to treat the physics of sur- 
face plasmons separately from the rest of the surface and 
bulk excitations, so that a field-theory of surface plas- 
mons in a curved space-time background may be con- 
sidered. For example, a nanohole in a thin metal mem- 
brane may be treated as a "wormhole" connecting two 
" fiat" surface plasmon worlds located on the opposite in- 
terfaces of the membrane [10]. On the other hand, near 
the plasmon resonance (which is defined by the condi- 
tion that em(w) — —Cd at the metal-dielectric interface, 
where em(a;) and Cd are the dielectric constants of metal 
and dielectric, respectively [9]) the surface plasmon ve- 
locity vanishes, so that the surface plasmon "stops" on 
the metal surface, and the surface charge and the nor- 
mal component of the electric field diverge. For a given 
frequency of light, the spatial boundary of the plasmon 
resonance ("the event horizon" of our toy model) may 
be defined at will using the geometry td{x,y) of the ab- 
sorbed layer of dielectric on the metal surface. Thus, 
the plasmon resonance becomes a natural candidate to 
emulate the event horizon of a black hole. As a result, 
toy two-dimensional surface plasmon black holes can be 
easily produced and studied [6,8]. 

In what follows I am going to explore the nonlinear op- 
tical behavior of the surface plasmon " black holes" based 
on the droplets of liquid on the gold or silver surfaces. I 
will show that such droplets exhibit strongly enhanced 
nonlinear optical behavior in the frequency range near 
the surface plasmon resonance of a metal-liquid inter- 
face due to considerable optical field enhancement near 
the toy event horizon (the edge of a droplet). This en- 
hancement may be responsible for the missing orders of 
magnitude of field enhancement in the surface enhanced 
Raman scattering effect [11]. In addition, experimental 
observation of the recently predicted [12] optical second 
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harmonic generation near the toy event horizon will be 
reported. Finally, the possibility of toy Hawking radi- 
ation observation will be discussed, and the expression 
for the effective Hawking temperature of a toy surface 
plasmon black hole will be derived. 

As a first step, let us consider the basic features of 
the surface plasmon toy black hole analogy. Let us con- 
sider in detail the dispersion law of a surface plasmon 
(SP), which propagates along the metal-dielectric inter- 
face. The SP field decays exponentially both inside the 
metal and the dielectric. Inside the dielectric the decay 
exponent is roughly equal to the SP wave vector. As 
a first step let us assume that both metal and dielectric 
completely fill the respective z <0 and z > half-spaces. 
In such a case the dispersion law can be written as [9] 

^^ = ^^$#., (1) 

where we will assume that = 1 — ujp/uj'^ accord- 
ing to the Drude model, and u>p is the plasma frequency 
of the metal. This dispersion law is shown in Fig. 1(b) 
for the cases of metal-vacuum and metal-dielectric in- 
terfaces. It starts as a "light line" in the respective di- 
electric at low frequencies and approaches asymptotically 
u! = Ljp/ (1 -I- Ed)^^^ at very large wave vectors. The latter 
frequency corresponds to the so-called surface plasmon 
resonance. Under the surface plasmon resonance condi- 
tions both phase and group velocity of the SPs is zero, 
and the surface charge and the normal component of the 
electric field diverge. Since at every wavevector the SP 
dispersion law is located to the right of the " light line" , 
the SPs of the plane metal-dielectric interface are decou- 
pled from the free-space photons due to the momentum 
conservation law. 

If a droplet of dielectric (Fig. 1(a)) is placed on the 
metal surface, the SP dispersion law will be a function 
of the local thickness of the droplet. Deep inside the 
droplet far from its edges the SP dispersion law will 
look similar to the case of a metal-dielectric interface, 
whereas near the edges (where the dielectric is thinner) 
it will approach the SP dispersion law for the metal- 
vacuum interface. As a result, for every frequency be- 
tween u>p/{l + e^)^/^ and ijjp/2^^'^ there will be a closed 
linear boundary inside the droplet for which the surface 
plasmon resonance conditions are satisfied. Let us show 
that such a droplet of dielectric on the metal interface be- 
haves as a "surface plasmon black hole" in the frequency 
range between Wp/(1 -I- ed)^/^ and LOp/2^/'^, and that the 
described boundary of the surface plasmon resonance be- 
haves as an "event horizon" of such a black hole. 

Let us consider a SP within this frequency range, which 
is trapped near its respective "event horizon", and which 
is trying to leave a large droplet of dielectric (see Fig. 2). 
The fact that the droplet is large means that ray op- 
tics may be used. Since the component of the SP mo- 
mentum parallel to the droplet boundary has to be con- 
served, such a SP will be totally internally reflected by 



the surface plasmon resonance boundary back inside the 
droplet at any non-zero angle of incidence. This is a sim- 
ple consequence of the fact that near the " event horizon" 
the effective refractive index of the droplet for surface 
plasmons is infinite (according to eq.(l), both phase and 
group velocity of surface plasmons is zero at surface plas- 
mon resonance). On the other hand, even if the angle of 
incidence is zero, it will take the SP infinite time to leave 
the resonance boundary. Thus, the droplet behaves as 
a black hole for surface plasmons, and the line near the 
droplet boundary where the surface plasmon resonance 
conditions are satisfied plays the role of the event horizon 
for surface plasmons. 

The toy black holes described above are extremely easy 
to make and observe [6]. In our experiments a small 
droplet of glycerin was placed on the gold film surface 
and further smeared over the surface using lens paper, 
so that a large number of glycerin microdroplets were 
formed on the surface (Fig. 3(a)). These microdroplets 
were illuminated with white light through the glass prism 
(Fig. 1(a)) in the so-called Kretschman geometry [9]. The 
Kretschman geometry allows for efficient SP excitation 
on the gold-vacuum interface due to phase matching be- 
tween the SPs and photons in the glass prism. As a re- 
sult. SPs were launched into the gold film area around the 
droplet. Photograph taken under a microscope of one of 
such microdroplets is shown in Fig. 3(b). The white rim 
of light near the edge of the droplet is clearly seen. It 
corresponds to the effective SP event horizon described 
above. Near this toy event horizon SPs are stopped or 
reflected back inside the droplet. In addition, a small 
portion of the SP field may be scattered out of the two- 
dimensional surface plasmon world into normal three- 
dimensional photons. These photons produced the image 
in Fig.3(b). We also conducted near- field optical mea- 
surements of the local surface plasmon field distribution 
around the droplet boundary Fig.3(c,d,e) using a sharp 
tapered optical fiber as a microscope tip. These measure- 
ments were performed similar to the measurements of 
surface plasmon scattering by individual surface defects 
described in [13]. The droplet was illuminated with 488 
nm laser light in the Kretschman geometry. The tip of 
the microscope was able to penetrate inside the glycerin 
droplet, and measure the local plasmon field distribution 
both inside and outside the droplet. Inside the droplet 
(in the right half of the images) the shear-force image (c) 
corresponds to the increase in viscous friction rather than 
the droplet topography. However, this image accurately 
represents the location of the droplet boundary, shown by 
the arrow in Fig. 3(e). The sharp and narrow local maxi- 
mum of the surface plasmon field just inside the droplet 
near its boundary is clearly visible in the near-field image 
Fig. 3(d) and its cross-section Fig. 3(e). 

Unfortunately, the described toy SP black hole model 
does not work outside the frequency range between 
ujp/il + SdY^"^ and ujpj'l}-!'^ . On the other hand, this is a 
common feature of every electromagnetic toy black hole 
model suggested so far. All such toy models necessarily 
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work only within a limited frequency range. In addition, 
the losses in metal and dielectric described by the so far 
neglected imaginary parts of their dielectric functions will 
put a stop to the singularities of the field somewhere very 
near the toy event horizon (however, such "good" metals 
as gold and silver have very small imaginary parts of their 
dielectric constants, and hence, very pronounced plas- 
mon resonances [9]). Notwithstanding these limitations, 
the ease of making and observing such toy SP black holes 
makes them a very promising research object. If we forget 
about the language of "black holes" and "event horizons" 
for a moment, the SP optics phenomenon represented in 
Fig. 3(b) remains a potentially very interesting effect in 
surface plasmon optics. Namely, this photo shows the 
existence of a two-dimensional SP analog of whispering 
gallery modes, which are well-known in the optics of light 
in droplets and other spherical dielectric particles. Whis- 
pering gallery modes in liquid microdroplets are known 
to substantially enhance nonlinear optical phenomena 
due to cavity quantum electrodynamic effects [14]. One 
may expect even higher enhancement of nonlinear opti- 
cal mixing in liquid droplets on the metal surfaces due 
to enhancement of surface electromagnetic field inher- 
ent to surface plasmon excitation, and in addition, due 
to accumulation of SP energy near the surface plasmon 
event horizons at the droplet boundaries. This strong en- 
hancement of nonlinear optical effects in liquid droplets 
may be very useful in chemical and biological sensing 
applications. It may also be responsible for the missing 
orders of magnitude of field enhancement in the surface 
enhanced Raman scattering (SERS) effect [11], since var- 
ious plasmon excitations are believed to play a major role 
in SERS. 

The current explanation of SERS is based on the 
combination of electromagnetic and " chemical" enhance- 
ments [11]. The current calculations of the electromag- 
netic field enhancement take into account "electrostatic" 
field enhancement at the apexes of various surface pro- 
trusions (the "lightning rod effect"), and the local field 
enhancement due to excitation of various localized sur- 
face plasmon modes in the crevices of the rough metal 
film. In addition, various weak and strong surface plas- 
mon localization effects are considered in combination 
with the consideration of a rough metal surface as a frac- 
tal object [15]. The "chemical" enhancement was pro- 
posed as an explanation for the quite a few missing or- 
ders of magnitude in theoretically calculated optical field 
enhancement, which follows from the magnitude of the 
experimentally measured SERS signals [11]. The "chem- 
ical" enhancement may happen if the molecular energy 
levels are affected by the proximity to the rough metal 
surface, and are drawn into resonance with the excitation 
field. 

SERS observations are usually conducted under not 
so well controlled conditions when the surface topogra- 
phy is not well-known and "dirty" (the target molecules 
are present on the surface in random locations). On the 
other hand, it is well known that even monolayer sur- 



face coverages considerably shift plasmon resonance [9] 
of the metal- vacuum interface. It would not be a big 
stretch to suggest that some areas of the " dirty" metal 
surface may contain compact areas covered with the mul- 
tiple layers of the target or solvent molecules (even if 
there are no droplets on the surface). Such compact ar- 
eas would affect surface plasmon propagation in a way, 
which is very similar to the effect of the droplets consid- 
ered above: effective surface plasmon event horizons may 
appear near the boundaries of such areas. Experimental 
data shown in Fig. 3 and the theoretical arguments above 
strongly indicate that the local field enhancement near 
these boundaries may be considerable. Even the data 
shown in Fig. 3(e) obtained with the limited optical reso- 
lution of the order of 100 nm [13] indicate at least 10-fold 
enhancement of the square of the local field intensity near 
the droplet boundary. In the physical picture of surface 
plasmon whispering gallery modes the local SERS sig- 
nal may be expected to grow by a factor of , where 
Q is the cavity quality factor [14]. For the surface plas- 
mons in the visible range Q may be estimated roughly as 
Q ~ L/2-kR, where L is the surface plasmon free prop- 
agation length and R is the droplet radius. Taking into 
account the typical theoretical value of L ^ AQjim [16] 
in the visible range, ^ 50 may be obtained for a 
R = 1/im droplet. Thus, such surface plasmon "whis- 
pering galleries" /" black holes" may provide considerable 
SERS enhancements on top of the local elec;tromagnetic 
enhancement due to other effects associated with the sur- 
face roughness. 

The effects of surface roughness on the properties of 
plasmons near the effective event horizon have not been 
considered so far. However, detailed analyses of these 
effects conducted recently for the case of real black holes 
[12] indicates that these effects will be especially strong in 
the nonlinear optics domain. For example, strong optical 
second harmonic generation has been predicted near the 
event horizon due to the effects of weak localization [12]. 

Wave propagation and localization phenomena in ran- 
dom media have been the topic of extensive studies dur- 
ing the last years [17]. One of the most striking exam- 
ples of such phenomena is the strong and narrow peak 
of diffuse second harmonic light emission observed in the 
direction normal to a randomly rough metal surface (see 
Fig. 4(a)). This peak is observed under the coherent illu- 
mination at any angle. This effect was initially predicted 
theoretically [18] and later observed in the experiment 
[19]. The enhanced second harmonic peak normal to the 
mean surface arises from the fact that a state of momen- 
tum k introduced into a weakly localized system will en- 
counter a significant amount of backscattering into states 
of momentum centered about -k. When these surface k 
and -k modes of frequency ui interact through an optical 
nonlinearity to generate 2a; radiative modes, the 2lo light 
has nonzero wave vector components only perpendicular 
to the mean surface. The angular width of the normal 
peak can be as small as a few degrees, and its amplitude 
far exceeds the diffuse omnidirectional second harmonic 
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background. 

Similar weak localization effects may be expected for 
the surface plasmon modes, which are trapped near the 
effective event horizon of our toy black hole when the 
metal surface exhibits moderate roughness. Similar to 
the case of planar rough surface [18,19], the momentum 
component parallel to the edge of the droplet must be 
conserved in nonlinear optical processes. This means that 
while the surface plasmon k and -k modes of frequency 

which interact through an optical nonlinearity, gener- 
ate 2<jo modes, the 2w plasmons could have nonzero wave 
vector component only perpendicular to the mean edge of 
the droplet (Fig. 4(b)). Thus, weak localization effects in 
the scattering of surface optical modes near the horizon 
should produce a pronounced peak in the angular dis- 
tribution of second harmonics of plasmons in the direc- 
tion normal to the horizon. In addition, because of such 
propagation direction (perpendicular to the edge of the 
horizon), these second harmonic plasmons have the best 
chances to escape the vicinity of the toy black hole. As 
a result, the relative intensity of the second harmonic ra- 
diation due to the weak localization effect near the event 
horizon will be much higher with respect to the diffuse 
SHG than in the case of planar rough surface. The dif- 
fuse SH plasmons would remain trapped near the horizon 
similar to the case of a real black hole [12]. 

Our experiments strongly indicate enhancement of 
SHG near the toy surface plasmon black holes. The 
droplets shown in Fig. 5(a) were illuminated by the 
weakly focused beam (illuminated spot diameter on the 
order of 50 /im) from a Ti:sapphire laser system con- 
sisting of an oscillator and a regenerative amplifier op- 
erating at 810 nm (repetition rate up to 250 kHz, 100-fs 
pulse duration, and up to 10 /iJ pulse energy), which was 
directed onto the sample surface in the Kretschman ge- 
ometry (Fig. 1(a)). Excitation power at the sample sur- 
face was kept below the ablation threshold of the gold 
film. The local SHG from the two droplets illuminated 
by the laser can be clearly seen in Fig. 5(b), which was ob- 
tained using a far-field microscope and a CCD camera. 
It should be noted, however, that while SH frequency 
plasmons experience the "event horizon" near the edges 
of the droplets, the refractive index of glycerine is not 
suflaciently large for the fundamental plasmons to see it. 
Thus, 810 nm plasmons may only be trapped into the 
whispering gallery modes near the edges of the droplets. 
On the other hand, the basic weak localization mecha- 
nism remains the same in the studied experimental situ- 
ation. 

Finally, let us examine the potentially most interesting 
nonlinear optical phenomenon, which may be exhibited 
by our toy black holes. The Hawking radiation [4] has 
been one of the most fascinating quantum physical phe- 
nomena discovered recently. However, nobody has ob- 
served it in the experiment so far, and its observation 
remains a great experimental challenge. In order to eval- 
uate the chances to see it in the experiments performed 
with surface plasmon toy black holes, we must first de- 



rive the expression for the effective Hawking tempera- 
ture. A derivation of an effective Hawking temperature 
of a toy dielectric black hole has been performed recently 
by Reznik [1] . Below we will closely follow his procedure. 

In order to emulate the effective space-time geometry 
considered in [1], let us consider a thin metal membrane 
with two "linear" droplets positioned symmetrically on 
both sides of the membrane (Fig.6(a)). The droplets 
thicknesses at a; = correspond to the surface plasmon 
resonance at the illumination frequency. The droplets ta- 
per off adiabatically in positive and negative x-directions 
on both sides of the membrane. In the symmetric mem- 
brane geometry the surface plasmon spectrum consists of 
two branches w_ and u)+, which exhibit positive and neg- 
ative dispersion, respectively, near the surface plasmon 
resonance [20]. Fig. 6(b) shows the dispersion curves of 
both branches for the cases of metal- vacuum interface far 
from the droplets, and for the locations near x = 0. The 
droplets represent an effective black hole for a;_ modes 
and an effective white hole for lo^ modes, similar to [1]. 
If the effective slowly varying local surface plasmon phase 
velocity c*{x) = ui/k{x) is introduced, the wave equation 
for surface plasmons may be written as 

Q2 Q2 Q2 

which corresponds to an effective metric 

ds^ = c*^dt^ - dx^ - dy'' (3) 

The event horizon corresponds to c* = at a; = 0. 

The two regions x > and a; < correspond to the two 
sides of a black hole, which are connected by an Einstein- 
Rosen bridge at a; = 0. The behavior of c*{x) near x = 
may be defined at will by choosing the corresponding ge- 
ometry of the droplet edge (if necessary, the droplet may 
be replaced by a similar shaped layer of solid dielectric). 
In order to adhere to the metric considered in [1] let us 
assume that c* = axe in the vicinity of x = 0. How- 
ever, we should remember that this linear behavior will 
be cut off somewhere near the effective horizon due to 
such effects as Landau damping, losses in the metal and 
the dielectric, etc. [6,8]. The resulting effective metric 
now looks like 

ds^ = a^x^c^df - dx^ - dy^ (4) 

Thus, this particular choice of the shape of the droplet 
edge gives rise to an effective Rindler geometry. Similar 
metric would describe the space-time geometry near the 
event horizon of a black hole with mass Mbh = c^/8ja, 
where 7 is the gravitation constant [1]. 

Following the standard derivation procedure [1,2,4,5], 
we must now introduce the analogue "Minkowski" coor- 
dinates as ?7 = — xe~"* and V = xe"*, which make the 
effective metric above conformally flat, and require that 
in the fiat "Minkowski" space the plasmon field will be 
in its vacuum state | Om >• However, the "Minkowski" 
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vacuum looks like a bath of thermal radiation with the 
temperature 

where Mp^ — {hc/jy^^ is the Planck mass and ks is 
the Boltzmann constant. For more details on the deriva- 
tion of this result one may address [1], where it is also 
demonstrated that cutting off c* near a; = and inclu- 
sion of the dispersion does not eliminate the Hawking 
radiation. 

We should point out that the final result for the effec- 
tive Hawking temperature docs not depend on the grav- 
itation constant 7. This is an encouraging fact since we 
are dealing with the toy black holes. In order to get 
a numerical estimate on the effective Hawking temper- 
ature we must assume some realistic value of a. Our 
approximation of a slow adiabatically changing c* inside 
the droplet means that c* does not change considerably 
on the scale of the local wavelength A* of surface plas- 
mons. Thus, a good top estimate for o: should look like 
a < 1/A ~ 1/Ao, where A is the plasmon wavelength far 
from the droplet, and Aq is the wavelength of light in 
vacuum. As a result, the effective Hawking temperature 
of our toy black hole is of the order of 

where u is the characteristic surface plasmon fre- 
quency. This value is quite close to the room temper- 
ature, which means that the droplet may probably con- 
form itself to be exactly in thermal equilibrium with its 
ambient. Thus, we come to a surprising conclusion that 
the shape of liquid droplets on metal surfaces may to 
some degree be determined by the Hawking radiation. 

However, we should accept this result with a degree 
of caution because of the recent claim that the electro- 
magnetic black hole analogues exhibit only the classical 
features of the black holes, while such quantum mechan- 
ical properties as Hawking radiation may not be repro- 
duced [21]. The basic reason for not having thermal ra- 
diation emitted is the energy conservation: there should 
be a source of thermal energy emitted to infinity. On 
the other hand, we must remember that the surface plas- 
mons are non-radiative modes, which arc bound to the 
interface. In addition, they do not propagate far alone 
the interface, being absorbed by the metal within a few 
tens of micrometers from the source. As a result, the 
surface plasmon Hawking radiation associated with the 
toy black hole would be observable only in the near field 
of the black hole, and energy conservation would not be 
violated. This situation would be rather similar to the 
recently discovered near-field thermal surface phonon- 
polariton emission in SiC [22], when the monochromatic 
emission with photon energies not represented in the far 
field zone may be observed in the near field of the SiC 
surface. In fact, this interesting effect is caused by the 



fact that c* ^ under the conditions of surface phonon- 
polariton resonance, similar to the toy Hawking radia- 
tion effect considered above. In addition, one may create 
non-equilibrium situations when energy is pumped into 
a toy black hole. Under such non-equilibrium conditions 
Hawking radiation may also be observable in the experi- 
ment. 

In conclusion, nonlinear optical behavior of recently 
introduced surface plasmon toy black hole models has 
been considered. Physical realization of these models in- 
volves droplets of liquid on a metal surface, which sup- 
ports propagation of surface plasmons. Such droplets arc 
shown to exhibit strongly enhanced nonlinear optical be- 
havior in the frequency range near the surface plasmon 
resonance of a metal-liquid interface. This enhancement 
may be responsible for the missing orders of magnitude of 
field enhancement in the surface enhanced Raman scat- 
tering effect. In addition, experimental observation of 
the recently predicted optical second harmonic genera- 
tion near the toy event horizon has been reported. Fi- 
nally, the possibility of toy Hawking radiation observa- 
tion has been discussed, and the expression for the effec- 
tive Hawking temperature of a toy surface plasmon black 
hole has been derived. 

This work has been supported in part by the NSF 
grants ECS-0210438 and ECS-0304046. 



[1] B. Reznik, Phys.Rev.D 62, 044044 (2000). 
[2] W.G. Unruh, Phys.Rev.D 51, 2827 (1995). 
[3] W. Sclileich and M.O. Scully, in New Trends in Atomic 
Physics, edited by G. Grynberg and R. Stora (North Hol- 
land, Amsterdam, 1984) p. 997. 
[4] S. Hawking, Nature 248, 30 (1974). 
[5] W.G. Unruh, Phys.Rev.D 14, 870 (1976). 
[6] I.I. Smolyaninov and C.C. Davis, gr-qc/0306089, submit- 
ted for publication in Phys. Rev. Lett. 
[7] L.D. Landau and E.M. Lifshitz, Field Theory (Pergamon, 

New York, 1984). 
[8] I.I. Smolyaninov, cond-mat/0309356, submitted for pub- 
lication in the New Journal of Physics. 
[9] H. Raether, Surface Plasmons, Springer Tracts in Mod- 
ern Physics Vol.111 (Springer, Berlin, 1988). 
[10] I.I. Smolyaninov, Phys.Rev.B 67, 165406 (2003). 
[11] A.G. Malshukov, Phys.Reports 194, 343 (1990). 
[12] I.I. Smolyaninov, Phys.Letters A 300, 375-380, (2002). 
[13] I.I. Smolyaninov, D.L. Mazzoni and C.C. Davis, 

Phys.Rev.Letters, 77, 3877 (1996). 
[14] A.J. 

Campillo, J.D. Eversole, and H.B. Lin, Phys.Rev.Lett. 
67, 437 (1991). 

[15] V.M. Shalaev, Nonlinear Optics of Random Media 

(Springer, Berlin, 2000). 
[16] P. Dawson, F. de Fornel, and J.-P. Goudonnet, 

Phys.Rev.Lett. 72, 2927 (1994). 



5 



[17] A. Ishimaru, Wave Propagation and Scattering in Ran- 
dom Media (IEEE, New York, 1997). 

[18] A.R. McGurn, T.A. Leskova, and V.M. Agranovich, 
Phys.Rev.B 44, 11441 (1991). 

[19] O.A. Aktsipetrov, V.N. Golovkina, O.I. Kapusta, T.A. 
Leskova, and N.N. Novikova, Phys.Lett.A 170, 231 
(1992). 

[20] K.L. Kliewer and R. Fuchs, Phys.Rev. 144, 495 (1966). 
[21] W.G. Unruh and R. Schutzhold, Phys.Rev. D 68, 024008 
(2003). 

[22] A.V. Shchogrov, K. .loulain, R. Carrninati, and J.-J. Gr- 
effct, Phys.Rcv.Letters 85, 1548 (2000). 

Figure captions. 

Fig.l (a) Experimental geometry of surface plasmon 
toy black hole observation, (b) Surface plasmon dis- 
persion law for the cases of metal-vacuum and metal- 
dielectric interfaces 

Fig. 2 A surface plasmon trapped inside a droplet near 
the effective "event horizon": The projection of surface 
plasmon momentum parallel to the droplet edge must 
be conserved. Due to effectively infinite refractive index 
near the droplet edge surface plasmons experience total 
internal reflection at any angle of incidence. 

Fig. 3 Far-field (a,b) and near-field (c,d) images of a toy 
surface plasmon black hole: (a) Droplet of glycerin on a 
gold film surface (illuminated from the top). The droplet 
diameter is approximately 15 micrometers, (b) The same 
droplet illuminated with white light in the Kretschman 
geometry, which provides efficient coupling of light to sur- 
face plasmons on the gold-vacuum interface (Fig. 1(a)). 
The white rim around the droplet boundary corresponds 
to the effective surface plasmon "event horizon", (c) and 

(d) show 10 X lO/itm^ topographical and near-field optical 
images of a similar droplet boundary (droplet is located 
in the right half of the images) illuminated with 488 nm 
laser light. Cross-sections of both images arc shown in 

(e) . Position of a droplet boundary is indicated by the 
arrow. 

Fig.4 (a) Strong enhancement of second harmonic 
emission from a randomly rough metal surface in the di- 
rection normal to the surface. A typical angular distri- 
bution of the diffuse second harmonic light is shown in 
the inset, (b) Second harmonic generation near the toy 
black hole horizon. Surface roughness is represented by 
a rough droplet edge. Second harmonic plasmons emit- 
ted perpendicular to the horizon have the best chances 
to escape the toy black hole. 

Fig. 5 (a) Microscopic image of the glycerin droplets 
under normal illumination, (b) SHG is seen from the 
droplets illuminated by 810 nm Ti:sapphire laser light in 
the Kretschman geometry shown in Fig. 1(a). 

Fig. 6 (a) A thin metal membrane with two "linear" 
droplets positioned symmetrically on both sides of the 
membrane. An effective event horizon is located at a; = 0. 
Such geometry emulates the effective space-time geome- 
try considered in [1]. (b) Dispersion laws of the lu^ and 
W-)- modes exhibit positive and negative dispersion, re- 
spectively, near the surface plasmon resonance. These 



branches arc shown for the cases of metal-vacuum in- 
terface far from the droplets, and for the locations near 
a; = 0. 
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